Stability of the Bi2Se3(lll) topological state: electron-phonon and -defect scattering 
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The electron dynamics of the topological surface state on Bi2Se3(lll) is investigated by 
temperature-dependent angle-resolved photoemission. The electron-phonon coupling strength is 
determined in a spectral region for which only intraband scattering involving the topological surface 
band is possible. The electron-phonon coupling constant is found to be A = 0.25(5), more than an 
order of magnitude higher than the corresponding value for intraband scattering in the noble metal 
surface states. The stability of the topological state with respect to surface irregularities was also 
tested by introducing a small concentration of surface defects via ion bombardment. It is found that, 
in contrast to the bulk states, the topological state can no longer be observed in the photoemission 
spectra and this cannot merely be attributed to surface defect-induced momentum broadening. 



PACS numbers: 73.20.At,71.70.Ej,79.60.-i 

Topological insulators are one of of the most intrigu- 
ing subjects of current condensed matter physics [T]-[3]. 
Despite of their insulating bulk, these materials support 
metallic edge and surface states with an unconventional 
spin texture [H [5] , electron dynamics [H [7] and stabil- 
ity. Exploiting these properties is the key to several ap- 
plications, e. g. in spintronics and quantum comput- 
ing. Moreover, several novel physical phenomena are pre- 
dicted in connection with the topological states |gHlQ]. 

The stable existence of a gap-closing surface state pTI - 
[T3] is a property derived from the bulk band structure of 
a topological insulator, rather than a mere coincidence. 
The topological state is also stable in a dynamical sense. 
A hallmark is the absence of back-scattering near non- 
magnetic defects. Edge states in the quantum spin Hall 
effect, a two-dimensional topological insulator, are com- 
pletely protected from (elastic) scattering [11 whereas 
the scattering phase space is strongly reduced for surface 
states on a three dimensional topological insulator [7|[T4], 
preventing localization by weak disorder. 

The stability of the topological state is essential for 
many phenomena and applications, however, only a few 
experimental studies have addressed this issue. Exper- 
imental measurements using angle resolved photoemis- 
sion spectroscopy (ARPES) and scanning tunnelling mi- 
croscopy (STM) have shown that the topological surface 
states are robust against a small number of adsorbates 
[151 m] and detectable at room temperature [15], but 
other results question their stability with respect to sur- 
face scattering processes [17 . Here we determine the 
electron-phonon {el — ph) coupling strength on the topo- 
logical insulator Bi2Se3(lll) [18l[T9]. In the absence of 
defects, el—ph scattering can be expected to be the domi- 
nant scattering mechanism at finite temperature and it is 
therefore of exceptional importance for any application. 



We concentrate on the spectral region in which only the 
topological state exists and thus only intraband scatter- 
ing is possible and we show that while the el — ph cou- 
pling constant A is of moderate size, it is surprisingly high 
compared to well-studied model systems such as the no- 
ble metal (111) surface states. We also demonstrate that 
scattering from a relatively small concentration of sur- 
face defects leads to a situation in which the topological 
surface state can no longer be observed by ARPES. 

ARPES experiments were performed on in situ- 
cleaved single crystals of Bi2Se3 at the ASTRID syn- 
chrotron radiation facility. Spectra were measured with 
an angular resolution of 0.13° and a combined energy 
resolution better than 15 meV. All data were taken at 
a photon energy of hu = 16 eV, except for the photon 
energy scan in Fig. 3 for which 14 < hu < 32 eV with a 
step size of Ahu = 0.1 eV. The biggest challenge in this 
type of study is the time-dependent band bending [151 [21] 
and the accompanying change of the surface electronic 
structure. This necessitates a fast measurement. Here 
the first ARPES measurements were typically taken 10 
minutes after cleaving the sample and the entire temper- 
ature dependence was concluded after another 5 minutes. 
Only the data shown in Fig. 3 were taken several hours 
after cleaving the sample. The STM measurements were 
performed at room temperature, also on samples cleaved 
in situ. More experimental details can be found in the 
supplementary online material (SOM) [22 . 

The electronic structure of the Bi2Se3 crystals is shown 
schematically in Fig. [ija). The topological state is situ- 
ated in the gap between the valence band (VB) and con- 
duction band (CB) and it shows the typical Dirac-cone 
dispersion. It is challenging to grow intrinsic samples of 
Bi2Se3 and it is often found, as in our case, that the ma- 
terial is degenerately n-doped with the Fermi level in the 
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bottom of the conduction band. This problem can be 
solved by counter-doping with Ca [23] but this has the 
disadvantage of introducing additional defects and it is 
therefore not a desirable strategy for the lifetime investi- 
gations performed here. 

Fortunately, however, the doping is not a problem for 
an investigation of the el—ph scattering of the topological 
state. For this state only the intraband scattering can be 
studied without involving any bulk states at all. To see 
this, consider the case of a decaying hole in the topolog- 
ical state in Fig. [TJa)-(c). The hole is filled by another 
electron with a phonon providing energy and momentum. 
The electrons available for filling thus have to originate 
from a very limited energy interval of hu above (below) 
the hole for the emission (absorption) of a phonon with 
this energy. As the maximum phonon energy in the ma- 
terial is quite small (22 meV [24]) this energy range is 
limited. As long as the initial state energy lies well within 
the bulk band gap only intraband scattering of the topo- 
logical state is allowed. Note that this is much different 
from electron-electron scattering where all electrons with 
a binding energy less than that of the hole, including the 
ones from the bulk conduction band, contribute to the 
lifetime. 

The decay of the hole can be expected to be lim- 
ited further by considerations involving the electron spin. 
The probability that an available electron fills the hole is 
greater if its spin orientation is similar to that of the hole 
state. The weights for these transitions are represented 
in a sketch of the F factor (see Ref. [25]) which is shown 
in Fig. [l] (c). Electrons which originate from the bright 
regions are more likely to fill the hole. 

Temperature-dependent ARPES measurements give 
straight-forward access to the el — ph coupling strength. 
Photoemission spectra taken between 90 K and 400 K are 
shown also in Fig. [ij A quantitative analysis of the data 
was performed by fitting constant energy slices through 
the spectra (momentum distribution curves) with two 
Lorentzian lines in order to obtain the position and width 
of the topological state over an energy range from 50 to 
100 meV above the Dirac point, well below the conduc- 
tion band minimum. The dispersion of the state was 
obtained from the peak position and fit with a polyno- 
mial to obtain a noise-free group velocity as a function 
of energy. This was combined with the width in order to 
determine the energy-dependent imaginary part of the 
self-energy (S^^. Not surprisingly, was found to be 
constant over the small energy range investigated so an 
average was taken to more precisely determine T^" . 

Figure |2] (a) shows the resulting as a function of 
temperature in sets corresponding to a number of dif- 
ferent samples and cleaves. The magnitude of varies 
greatly depending on the sample and cleave quality which 
is a result of different and inhomogeneous defect concen- 
trations. What is also apparent is that has a con- 
sistent dependence on temperature. To see the temper- 
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FIG. 1: (color online), (a) Band structure of Bi2Se3 showing 
the topological surface state. A photohole of energy E can 
be filled by electrons which absorb or emit a phonon and 
originate from an energy window of zLhuu where hu is the 
maximum phonon energy in the material, (b) Possible el — ph 
scattering processes to fill a hole, (c) Probability that an 
electron with a certain k\\ will fill the photohole (bright areas 
having the highest probability), (d)- (f) ARPES spectra taken 
at different temperatures. 
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FIG. 2: (color online), (a) Imaginary part of the self energy 
(S^^) as a function temperature where each set of two points 
represents a different sample or cleave. The estimated un- 
certainty for each data point is =b 3 meV. (b) Data from (a) 
where the data sets with open circles have been shifted in en- 
ergy such that a line fit through all the data minimizes x^- 
The dotted line show a calculation for a full Debye model (see 
text). 
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ature dependence more clearly, the data sets displayed 
with open circles in Fig. [2] (a) have been rigidly offset 
in energy such that a line fit through all the data points 
minimizes (see Fig. [2](b)). 

The el—ph coupling strength can be expressed in terms 
of the el—ph coupling constant A. In the present context, 
A is a spectroscopic parameter defined via the Eliashberg 
function of the system rather than through the increase 
of the effective mass at the Fermi energy. At high tem- 
peratures, independent of phonon spectrum, A can be ex- 
tracted from the temperature-dependent using the re- 
lation A = {TikB)~^dT." /dT. While this is formally valid 
only at high temperatures, it is expected to hold rea- 
sonably well for temperatures comparable to the Debye 
temperature [26 and thus also in our case (Od = 182 K 
[27j). A line fit to the data results in A = 0.25±0.05. The 
actual value might be slightly smaller due to the effect of 
the time-dependent band bending [22 . 

To confirm the justification of the simple linear fit, we 
also calculate T," in the full Debye model using 6^ = 
182 K and A = 0.26 [28]. The result is shown as a dotted 
line in Fig. [2](b) along with the linear fit. The difference 
is very small and the change of A to get a similar slope 
at higher temperatures is well within the uncertainties. 
Note that the real difference is probably even smaller be- 
cause the Debye temperature at the surface is often only 
^70% of the bulk value [29^. In order to fit the Debye 
model to the absolute value of T^" reported in Fig. [2](b), 
we have to assume a temperature-independent contribu- 
tion to the self-energy of 20 meV, accounting for both 
electron-electron and electron-defect scattering. Thus, 
the el — ph coupling contribution is rather significant, 
amounting to half of the total self-energy at room tem- 
perature. Indeed, its relative significance can be expected 
to be higher still for intrinsic Bi2Se3 where the electron- 
electron contribution to the self energy should be much 
smaller than in our strongly n-doped sample. 

At first glance, a value of A = 0.25 appears to imply 
a moderate el — ph scattering for the topological state. 
The value is higher than for a good conductor such as the 
noble metals (A ~ 0.1) but much smaller than for a BCS 
superconductor (A ^ 1) [30]. For the particular system of 
a topological surface state, however, the value is surpris- 
ingly high because the phase-space for the scattering is 
very small. Only intraband scattering is allowed and the 
number of available states is further restricted by the re- 
quirement for spin conservation outlined above. Rather 
than to compare to the absolute A value of other systems, 
one should only consider the intraband contribution for 
a two-dimensional state. This is well understood for the 
(111) surface states of the noble metals for which the in- 
traband contribution to A is only ^ 0.02 [20 . Our results 
thus imply a very significant el — ph scattering strength 
for the topological state. This is a rather unexpected 
finding and in strong contrast to an earlier analysis of 
the topological state's self-energy where the seemingly 



reasonable assumption of a vanishing intraband scatter- 
ing was made [37. 

We briefly address the related issue of scattering by 
static defects. To this end, defects were created by using 
a 15 second cycle of sputtering with 0.7 keV Ne+ at a 
pressure of 3 x 10~^ mbar. The introduction of a small 
concentration of surface defects due to the sputtering is 
expected to strongly influence the two-dimensional (2D) 
states associated with the surface while they should only 
have a minor effect on the three-dimensional (3D) bulk 
states. 

Figure [3| summarizes the photoemission and STM re- 
sults for the pristine and sputtered surfaces. 2D and 3D 
states are best distinguished by the absence of kz dis- 
persion of the former. Therefore results from multiple 
photon energies are shown and the kz dependence of the 
band structure is calculated assuming a free-electron fi- 
nal state model [181 EI] • There are three features that do 
not disperse with respect to /c^, thus confirming their 2D 
character: the topological state, an "M" -shaped feature 
at a binding energy of ^ 0.75 eV and a surface-localized 
2D electron gas formed in the quantum well between vac- 
uum and the downwards-bent bulk conduction band [21] . 
The latter feature is only observed some time after the 
cleave and thus neither present in the data of Fig. [l] 
nor in that used to extract the self-energies in Fig. [2j 
A comparison of the data before and after sputtering 
shows a removal of the 2D features whereas the 3D fea- 
tures are merely broadened, most probably due to the 
non-conservation of /cy. 

The disappearance of both the 2DEG and the "M"- 
state can easily be attributed to the perturbation of the 
Bi2Se3 surface. What is much more difficult to account 
for is the disappearance of the topological surface state 
which, due to its topological nature, is expected to be ro- 
bust against surface defects such as those created by sput- 
tering. The surface defects can act as scattering centers 
and reduce the mean free path X^q of the surface state 
electrons, leading to an uncertainty in the crystal momen- 
tum on a scale of A/cy = l/A^gp [32 . If the surface defect 
density is high enough, the disappearance of the topolog- 
ical state could be explained by momentum broadening. 
We estimate the average separation between surface de- 
fects from STM images taken prior to (after) sputtering 
to be 425 A (85 A) (Fig. [3|. If every defect acts as a scat- 
tering center we thus obtain a possible crystal momentum 
broadening of Ak\\ = 0.002 A"^ and Ak\\ = 0.012 A-^ 
for the pristine and sputtered surfaces, respectively. A 
momentum broadening of Ak\\ = 0.012 A~^ would be far 
too small to explain the extinction of the topological sur- 
face states. Hence, it appears that the topological states 
are strongly perturbed, even by a small concentration of 
defects, consistent with recent transport measurements 
[TT] . Note that this observation is not necessarily incon- 
sistent with previous findings of stability after a long time 
in vacuum [31 or even exposure to air ^33^ if one assumes 
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FIG. 3: (color online). The effect of defect creation on the surface. Differences between (a) a pristine surface and (b) a sputtered 
surface. Each panel shows a three-dimensional representation of the electronic structure obtained by taking data as a function 
of photon energy. Images corresponding to the center of the bulk Brillouin zone are shown separately. The STM images of the 
surface were measured using a voltage and tunneling current of 1.84 V and 0.18 nA in (a) and 1.68 V and 0.45 nA in (b). 



the surface of Bi2Se3 to be chemically inert. Note also 
that the failure to observe the state in ARPES does not 
imply that it does not exist anymore either. Further re- 
search, including more local tools such as spectroscopic 
STM are needed to investigate this issue. 

In conclusion, we find that the el —ph coupling for the 
topological state of Bi2Se3 is surprisingly strong in view 
of the fact that only intraband scattering can contribute 
to the decay of the photohole. On the other hand, the 
state still dominates the spectral function at room tem- 
perature such that the el —ph coupling is not expected to 
preclude room temperature applications. The states' ap- 
parent sensitivity towards a moderate amount of defects 
is a more significant consideration for practical applica- 
tions. Despite the relatively inert surfaces of the topolog- 
ical insulators, it appears that a useful practical device 
will exploit states located at an interface between a topo- 
logical and a trivial insulator rather than on the interface 
between a topological insulator and vacuum. Such a con- 
figuration would have the added advantages that it would 
stabilize the fragile band bending near the interface and 
permit a gating of the near surface electronic structure. 
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